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Abstract 

The effective interaction between a quark and an anti-quark as obtained previously 
with by the method of iterated resolvents is replaced by the ||-model and applied 
to flavor off-diagonal mesons including the vr^. The only free parameters are the 
canonical ones, the coupling constant and the masses of the quarks. 



The light-cone approach [1] to the bound-state problem in gauge theory [2] 
aims at solving Hlc\'^) = M^l^f). If one disregards possible zero modes [4] 
and works in the light-cone gauge, the (light-cone) Hamiltonian H^c is given 
as a Fock-space operator in [2], or in the compendium [3]. Its eigenvalues 
are the invariant mass-squares of physical particles associated with the 
eigenstates In general, they are superpositions of all possible Fock states 
with its many-particle configurations. For a meson, for example, holds 
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If all wave functions like or "^gg are available, one can analyze hadronic 
structure in terms of quarks and gluons [2]. - Hiller [5], or Tang et al. [6], 
attacks this problem by DLCQ. Alternatively, one addresses to reduce the 
many-body problem behind a field theory to an effective one-body problem. 
The derivation of the effective interaction becomes then the key issue. By def- 
inition, an effective Hamiltonian acts only in the lowest sector of the theory 
(here: in the Fock space of one quark and one anti-quark) and has the same 
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eigenvalue spectrum as the full problem. I have derived such an effective in- 
teraction hj the method of iterated resolvents [7], that is by systematically 
expressing the higher Fock-space wave functions as functionals of the lower 
ones. In doing so the Fock-space is not truncated and all Lagrangian symme- 
tries are preserved. The projections of the eigenstates onto the higher Fock 
spaces can be retrieved systematically from the gg-projection, with explicit 
formulas given in [8]. I have derived [9] the effective interaction also with the 
method of Hamiltonian flow equations [10]. 



1 The eigenvalue equation in the qq-spa.ce 



For flavor off-diagonal mesons (to mesons with a different flavor for quark and 
anti-quark), the effective one-body equation [7] is shockingly simple: 
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Here, is the eigenvalue of the invariant-mass squared. The associated eigen- 
function ip = is the probability amplitude {x, k±; Xi, X2\'4') for finding a 
quark with momentum fraction x, transversal momentum k± and helicity Ai, 
and correspondingly the anti-quark with 1 — x, —k± and A2. The mi and m2 
are (effective) quark masses and a is the (effective) coupling constant. The 
mean Feynman-momentum transfer of the quarks is denoted by Q^, 
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the spinor factor 5* = S{x, k±; x', k[) by 

(Ai, X2\S\X[, X'2) = [u{ki, Xi)ru{k[, a;)] [v{k',, X',h,v{k2, X2)] . (3) 



The regulator function R{x', k'j_; A) restricts the range of integration as func- 
tion of some mass scale A. I happen to choose here a soft cut-off (see below), 
in contrast to the previous sharp cut-off [11,12]. Note that Eq.(l) is a fully 
relativistic equation. 

The full QCD-Hamiltonian in 3+1 dimensions must be regulated from the 
outset. There are arbitrarily many possible ways to do that, see f.e. Hiller [5]. 
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One of the few practical ways is vertex regularization [2,7], where every Hamil- 
tonian matrix element, particularly those of the vertex interaction (the Dirac 
interaction proper), is multiplied with a convergence-enforcing momentum- 
dependent function, which is kind of a form factor [2]. The precise form of this 
function is unimportant here, provided it prevents that the scattered particles 
go too much off-shell as function of a cut-off scale (A). In the limit A — > oo 
the full unregulated theory is restored. The effective quark masses mi and 
m2 and the effective coupling constant a depend, in general, on A. Explicit 
expressions for a sharp cut-off are available [7]. In the spirit of renormalization 
theory they are renormalization constants, subject to be determined by exper- 
iment, and hence-forward will be denoted by mi, m2, and a, respectively. In 
next-to-lowest order of approximation the coupling constant becomes a func- 
tion of the momentum transfer, a — > A), with the explicit expression 
given in [7]. 

Why is Eq.(l) so shocking? 

In the first place, I am baffled by the fact that the effective interaction for 
QCD and QED differs only by the color factor (n^ - l)/(2nc) = 4/3. This 
is in raging conflict with all what 1 have learned over the past 20 years on 
vacuum structure, condensates, soft pion theorems, conflnement, and the like. 
In fact, I knew the simple color factor already at the flrst light-cone meeting in 
1991. I have hidden it quite well in Ref.[ll], where I had obtained Eq.(l) by a 
light-cone adapted Tamm-Dancoff approach, which is an educated guess at the 
most. Like everybody else, I could not believe that Eq.(l) is applicable to the 
pion, for example, although it really should, since no assumption was made on 
the size of the quark masses. In consequence, I have spent in vain many painful 
years to accumulate counter evidence by a more rigorous formalism, first by 
investigating the zero modes [4], then by developing the method of iterated 
resolvents [7]. For some time around the 1998 meeting in St. Petersburg, I 
believed that the salvation was in the functional behaviour of the nmning 
coupling constant a[Q;h.). But without screwing the parameters beyond all 
reason, I could not get the pion quantitatively (unpublished), the pion that 
mystery particle of QCD. Ultimately I am thrown back to Eq.(l), the point 
of departure. 

But nowadays I am in a better shape. I have obtained Eq.(l) by two dif- 
ferent but well founded methods: by the iterated resolvents [7] and by the 
Hamiltonian fiow equations [10]. 

It might be to early for solving Eq.(l) numerically in full glory like in Ref.[12]. 
Rather should I try to dismantle the equation of all irrelevant details (or what 
I think they are) and work as much as possible analytically. With other words, 
I should develop a simple model. 
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2 The tJ,-model and its application to the pion 



In light-cone parametrization, the quarks are at relative rest when k± — and 
X — X = mi /{mi + 7712). For very small deviations from these equilibrium 
values the spinor matrix is proportional to the unit matrix, with 

(Ai, X2\S\X[X'2) ~ 4mim2 Sx,,x[ h^,x'^. (4) 



see Compendium. For very large deviations, particularly for k'^ » kl, holds 
Q^-k'^, and (Ti 1^1 Ti) =^ 2^1'- (5) 



Both extremes are combined in the "tJ,-model": 

(ti l-^l ti) _ (IT 1^1 it) 4mim2 



g2 - Q2 +2- (6) 



It interpolates between two extremes: For small momentum transfer, the '2' 
is unimportant and the dominantly Coulomb aspects of the first term prevail. 
For large momentum transfers the Coulomb aspects are unimportant and the 
hyperfine interaction is dominant. But the model over-emphasizes many as- 
pects: It neglects the momentum dependence of the Dirac spinors and thus 
the spin-orbit interaction; it also neglects the momentum dependence of the 
spin-spin interaction, keeping only the '2' as its residue. This '2' creates havoc: 
Its Fourier transform is a Dirac-delta function with all its consequences in a 
bound-state equation. 

Here is an interesting point: We are all familiar with the field theoretic diver- 
gences residing in the effective masses and the effective coupling constant. We 
are not used to "divergences" residing in a finite number '2'. They must also 
be regulated and renormalized. 

I regulate thus the kernel by 

5 / 4mim2 \ 4mim2 , 2/^^ 



I denote the soft cut-off by (in contrast to a sharp cut-off A). The Coulomb 
interaction in first term needs no regularization. In consequence I replace 
Eq.(l) by 
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where '4'{x, kj_) = {x, k±; t, i \i>)- On top of the canonical parameters m and 
a, the eigenvalues depend now on a regularization scale /i. 

Since /i is an unphysical parameter, I must remove its impact. I do that in 
the spirit we had treated the pairing interaction in nuclei [13]. It also has a 
Dirac-delta function. As emphasized repeatedly by the late and un-forgotten 
Strutinky [13], one must adjust the (pairing) coupling constant and the (pair- 
ing) cut-off such, that the binding energy and the pairing gap would not de- 
pend on either of them. I proceed therefore by looking for a cut-off dependent 
coupling function such, that the calculated mass of the pion agrees with 
the empirical value. The so obtained function is considered universal. The 
actual value of ^ must be fixed by a second requirement. I emphasize that this 
procedure is not a renormalization group treatment in the modern sense; but 
it is a first step in the right direction. 

For carrying out this programme in practice, I need an efficient tool for solv- 
ing Eq.(8). Such one has been developed recently [14]. I outline in short the 
procedure for the special case mi — m2 = m. 1 change integration variables 
from X to kz by the Sawicki transformation and substitute 

x{kz) = ^ + I =, k±) = , ^ (9) 

2 2^m^ + kl + kl 1 - x) 



The variables {k^, k±) are collected in a 3-vector p, with all further details in 
the Compendium [3]. Eq.(8) becomes then 



M'(l)(p) = [Am^ + 4p2J (10) 

" 3^ I"^'^' {{p-pV ^ m^^ + {p-p'y) 

For the present purpose it suffices to restrict to spherically symmetric s— states 
and to apply Gaussian quadratures with 16 points. On an alpha work station 
it takes a couple of micro-seconds to solve for the spectrum of a particular 
case. Fixing the up and the down mass to m„ = rud — m (the codes work with 
dimensionless units) 

m = 1.16 K = 406 MeV, k = 350 MeV, (11) 
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Fig. 1. Coupling constant a versus the Fig. 2. Spectrum of vr^: The first three 
regularization scale /i/k. singlet-s states in MeV versus /x/k. 



the calculated pion mass can be made to agree within eight digits with the 
experimental value [15]. For every value of /i one obtains thus an which is 
displayed in Fig. 1. For every value of /i (thus a^), I get a whole spectrum, as 
plotted in Fig. 2. The lowest state corresponds to the the vr"*" and stays nailed 
fixed to the empirical value. Much to my surprise, the second and the third 
(as well as the higher ones) expose an extremum. At the extremum, I fulfill 
Strutinsky's requirement 



= 0, with /io ^ 3.8 K, = 0.6904. (12) 



The regularization scale /i determines itself from the solution! 

Here comes the next problem: How should I determine the quark mass? In 
Eq.(ll), I selected a particular value of m; I could have taken any other and 
would have gotten ii{rn) and a^{rn). I obviously need a second empirical da- 
tum. Which one? I could choose the first excited state of tt"*"; but that is not 
known with sufficient accuracy for my purposes and, moreover, there are un- 
certainties of interpretation, see Morningstar [17] and Llanes- Estrada [18]. I 
could take the mass of the p^; but the present ||-model is not particularly 
suited for vector mesons. I remain thus with the pions size. Knowing the wave- 
function, see below, I could calculate the form factor [2] and thus the exact 
root-mean-square radius. But I replace this non-trivial calculation by the the 
following very drastic but simpler construction: I parametrize the numerical 
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Fig. 3. The pion wave function ^{p) 
is plotted versus p/(1.552k) in an ar- 
bitrary normalization. The filled circles 
indicate the numerical results, the open 
circles the fit function. The diamonds 
denote the pure Coulomb solution. 

results by a fit function with one free 
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Fig. 4. The relativistic potential energy 
W{r) for the pion is plotted in GeV ver- 
sus the radius in fm (solid line). The 
pion mass is indicated by the dashed, 
and the ionization threshold by the 
dashed-dotted line. See Eq.(16). 

parameter (pa) to be, 

with Pa = 1.471 K. (13) 



The quality of this fit can be judged by Fig. 3. 1 Fourier-transform this function 
to configuration space, giving tp{f) ~ e~^°'', and calculate trivially 



. x/3 

< >|= — = 0.663fm. (14) 

^ Pa 



Finally, I vary m until I get about agreement with < >2= 0.657 fm, the 
experimental value [15]. This exhausts all freedom in choosing the parameters 
of Eq.(8). Of course, the drastic assumption behind Eq.(14) has to replaced in 
the future by the correct prescription [20]. Actually, I had started the proce- 
dure with m = 1.00 K. Eq.(ll) quotes the result after iteration. After inserting 
the inverse Sawicki transformation Eq.(9) into the fit Eq.(13), I obtain 



x(l 



X 



rri^ {2x - if + k'l 



(15) 
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This completes my goal: I have a pion with the correct mass and size, and I 
have an analytic expression for its light-cone wave function. 

For the first time in my life I see in Eq.(15) a genuine light-cone wave function, 
which is directly related to the QCD-Lagrangian (although only after quite a 
few approximations and simplifications). It could be used thus as a baseline for 
calculating the higher Fock-space amplitudes, as explained in [8]. It differs from 
the literature [19] by the preceding factor. 1 like to emphasize that regulated 
delta-interaction, the Yukawa potential in Eq.(8) at the scale /iq, acts like 
a genuine Dirac-delta function: it pulls down essentially only one state, the 
pion, in analogy with the pairing model [13]. The others are left more or less 
unperturbed at their Bohr values, see Fig. 2. 1 like to emphasize further that 
1 am not in conflict with the usual perturbative inclusion of the hyperflne 
interaction. For sufficiently large values of /i, say for /i » 10, see Fig. 2, the 
coupling constant decreases strongly and the spectrum becomes more and 
more the familiar Bohr spectrum with a small hyperflne shift, such that it can 
be calculated perturbatively, indeed. In a future and more rigorous solution 
of the full equation 1 expect that the excited states can be disentangled into 
almost degenerate singlets and triplets, which in turn can be interpreted as 
an excitation of the pseudo-scalar pion, or the ground state of a vector meson, 
respectively. In any case, the flrst excited state of the simple model correlates 
very well with the mass of the p"*" and the other vector mesons, see below. 



3 Extensions and conclusions 

I can apply the same procedure also to the other mesons. Once I have the 
up and down mass, I determine the strange, charm and bottom quark mass 
by reproducing the masses of the K~ and respectively. This gives 
rus — 508, rric — 1666 and rri}, — 5054 MeV. All these numbers as well as 
those in Tables 1 and 2 are rounded for convenience; the calculations or the 
experiments are available with greater precision. This flxes all parameters 
which there possibly are. 1 generate with them all off-diagonal pseudo-scalar 
mesons and compile their mass in Table 1. The corresponding wave functions 
are also available but not shown. In view of the simplicity of the model, the 
agreement with the empirical values in Table 2 is remarkable. The mass of the 
flrst excited states in Table 1 correlates astoundingly well with the associated 
vector meson in Table 2. 

After so many years, I am forced to conclude that the pion looks quite a bit 
different than told in the literature. Although the above picture is kind of 
wood-carved, the previous work with Kalloniatis and Pinsky lets appear it 
unlikely that the zero modes, once included, change the picture drastically. I 
have found no evidence that the vacuum condensates are important. I conclude 
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Table 1 

The calculated mass eigenvalues in MeV. 
Those for singlet-ls states are given in 
the lower, those for singlet-2s states in 
the upper triangle. 
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871 


2030 


5418 
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494 




2124 


5510 
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1865 


1865 


1929 




6580 
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5279 


5279 


5338 


6114 





Table 2 

The empirical masses of the fiavor-off- 
diagonal physical mesons in MeV. The 
vector mesons are given in the upper, the 
scalar mesons in the lower triangle. 

u d 's c h 

u 768 892 2007 5325 

d 140 896 2010 5325 

s 494 498 2110 — 

c 1865 1869 1969 — 

b 5278 5279 5375 — 



that the pion is describable by a QCD: The very large coupling constant in 
conjunction with a very strong hyperfine (spin-spin-) interaction makes it a 
ultra strongly bounded system of constituent quarks. More then 80% of the 
constituent quark mass is eaten up by binding effects. No other physical system 
has such a property. But this pion is not the one I was dreaming of when this 
work began. 

The theory has three parameters. Two are determined by the pions mass and 
size, and one determines itself from the theory. In the lack of further empirical 
data, I have no independent check for the above strong statements, except 
that the first excited state should roughly be degenerate with the mass of the 
p-meson. At present I am working on a check [20] whether Eq.(15) is consistent 
with Ashery's experiment [21]. 

The question of confinement is closely related to a the potential, the strong 
potential K(r), but in the present approach the relation is subtle. One of the 
great advantages of light-conc quantization is the additivity of the free part 
and the interaction. Because of the dimensions of invariant mass squares the 
relation of potential and kinetic energy is somewhat hidden. What I can do, 
however is to Fourier transform Eq.(8) and to interpret the operators p and q as 
the classical analogue of momentum and position. Instead of an Hamiltonian 
H{q,p) I then have an invariant mass-squared M'^{q,p), but for the argument 
that does not matter. I thus have 

M(r,p)^mJl + + — Vsir) = ^(l + -^e-^A .(16) 

V msTTir rus or \ mstrir J 

The expansion into rest mass -|- kinetic energy -|- potential energy, 

2 

M{r,p) c^ms + ^ + Vs{r), (17) 
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is justified however only for <C mgrnr, and Vs{r) ^ m^, wliicli can not 
be satisfied for all r. I therefore propose to introduce a relativisiic poteniial 
energy by W{r) = M{r,p = 0) as the closest relativistic analogue to a non- 
relativistic potential Vs{r), which we all have in the back of our mind. This 
function is plotted in Fig. 4. It vanishes at ro = 0.5194 fm and the classical 
turning point is = 0.5300 fm. Since the sum of the quark masses is 812 MeV, 
the ionization threshold occurs at 542 MeV: In order to liberate the quarks in 
the pion, one has to invest more than three times the rest mass of a pion; I 
call this practical confinement. 



I contrast these results to Lattice Gauge Calculations. It is not generally known 
that LGC's have considerable uncertainty to extrapolate their results down to 
such light mesons as a pion; it exceeds their grid sizes. It is also not generally 
known that Lattice Gauge Calculations get always strict and linear confine- 
ment even for QED, where we positively know the ionization threshold. The 
'breaking of the string', or in a more physical language, the ionization thresh- 
old is one of the hot topics at the lattice conferences [22] . Moreover, in order 
to get the size of the pion, thus the form factor, requires another generation 
of computers and physicists to run them. 

I contrast these results also to phenomenological approaches. They usually 
do not address to get the pion, and for the heavy mesons, where they are 
so successful, a good phenomenological model has quite many parameters, 
in any case more that the above canonical ones. A detailed comparison and 
systematic discussion of the bulky literature can however be postponed, until 
we are ready to solve the full Eq.(l) and do not restrict ourselves to the 
caricature of the QCD- interaction, the ti-model. 

I contrast these results, finally, to the Nambu-Jona-Lasio-based models which 
are so successful in accounting for the isospin. I do not even think of quoting 
from the huge body of literature but mention in passing that the NJL-models 
are not renormalizable, have no relation to QCD, and deal mostly with the 
very light mesons. They just about can handle strangeness and break down 
for the heavy flavors. 

In conclusion I may state that not a single model can describe quantitatively 
all mesons from the tt up to the T from a common point of view, from QCD. 
By solving the tJ,-model I provide at least some evidence that Eq.(l) can 
eventually give that, perhaps sometimes in the future. 

Acknowledgement. I thank Susanne Bielefeld for her unselfish help during 
the past three years. 
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